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Abstract. Consider Dyson's Hermitian Brownian motion model after a finite 
time S, where the process is started at N equidistant points on the real line. 
These N points after time S form a determinantal process and has a limit as 
TV — > oo. This limting determinantal proceess has the interesting feature that 
it shows number variance saturation. The variance of the number of particles 
in an interval converges to a limiting value as the length of the interval goes to 
infinity. Number variance saturation is also seen for example in the zeros of the 
Riemann C-function, ISJ, 0. The process can also be constructed using non- 
intersecting paths and we consider several variants of this construction. One 
construction leads to a model which shows a transition from a non-universal 
behaviour with number variance saturation to a universal sine-kernel behaviour 
as we go up the line. 



The Bohigas-Gianonni-Schmidt conjecture, Uj|, says that the spectrum {Ej}j>i, 
Ej — ► oo as j — > oo , of a quantum system whose classical dynamics is fully chaotic, 
has random matrix statistics in the large energy limit, Ej — > oo. For finite Ej 
the specrum has non-universal features depending on the particular system, but as 
we go higher up in the spectrum the statistical properties become more and more 
like those from the universal point processes obtained from random matrix theory. 
The zeros of Riemann's C-function, {-Ej}j>i, < E\ < E% < . . . (assuming the 
Riemann hypothesis) show a similar behaviour and has been popular as a model 
system in quantum chaos, [3], since there are many analogies. The number of zeros 
< E, denoted by Af(E), is approximately M- log M- — If we unfold the zeros 
by letting Xj = Af(Ej), so that the mean spacing becomes 1, it is conjectured 
by Montgomery, [I^], and tested numerically by Odlyzko, [201 , that the statistics 
of the Xjis as j — > oo, is like the statistics of a determinantal point process with 
correlation functions 



m > 1. Weak forms of this have been proved in ^2], |21j . 

This determinantal point process is obtained as a scaling limit of GUE or the 
classical compact groups, e.g. U(n), and also as the universal scaling limit of many 
other hermitian random matrix ensembles. If we count the number of particles 
(eigenvalues) in an interval of length L in a process on M with correlation functions 
(II. 1|) we get a random variable with variance, the number variance, ~ 4ylogL as 
L — > oo (an exact formula for finite L is given in (|2.33|) below). A feature in the 
quantum chaos model and for the Riemann zeros is number variance saturation, |2j , 
\2'2\ . If we consider the Riemann zeros and intervals of length L at height E, where 



1. Introduction 



(1.1) 




2 



K. JOHANSSON 



L -C E, and compute the variance by considering many disjoint intervals of length 
L, the dependence on L is such that for small L < log it behaves like ^-logL 
but as L grows it saturates, actually oscillates around an average value which is 
approximately -\ log(log see the work of Berry, [2|, for interesting precise pre- 
dictions. Hence the sine kernel determinantal point process is only a good model in 
a restricted range which becomes longer as we go up the line. The question that we 
address in this paper is whether it is possible to construct a determinantal process 
which shows number variance saturation? Can we construct a determinantal pro- 
cess on [0,oo) which shows a transition from a non-universal regime to a universal 
regime described by as we go further and further away from the origin? It is 
not possible to get number variance saturation with a translation invariant kernel, 
like in i|l.lfl , since the sine kernel is the kernel with the slowest growth of the num- 
ber variance among all translation invariant kernels which define a determinantal 
point process, |2"3] . 

In this paper we will construct models having these properties by suitable scaling 
limits of determinantal process defined using non- intersecting Brownian motions. 
These models will not be translation invariant. We can restore translation invari- 
ance by averaging, but then we will no longer have a determinantal point process. 
One of the kernels obtained is given approximately by 

,j ^, sin7r(izi — y) ^ dcosir{x + y) + (y — x) sin7r(y + x) 

tt(x — y) tt(cP + (y — x) 2 ) 

where d > is a parameter, see (|2.23|) below. (The exact kernel has corrections 
of order exp(— Cd)). This model will have a number variance with saturation level 

~£log(27Td). 

There are other connections between .L-functions and random matrix theory. 
Katz and Sarnak, |15j , study low- lying zeros of families of L- functions and connect 
their statistical behaviour with that obtained for the eigenvalues of random matrices 
from the compact classical groups with respect to Haar measure. This leads to a 
classification of the L-functions into different symmetry classes. Three different laws 
for the distribution of the lowest zero are obtained. We will see below that these 
three laws can also be obtained from the non-intersecting Brownian motions by 
choosing different boundary conditions, see also ^DJ- Another recent development 
is the study of characteristic polynomials of matrices from the classical groups which 
have been used to model L-functions, and led to interesting conjectures for their 
moments, see [T^]. See also jS] for a discussion of linear statistics of zeros. 

In this paper we will have nothing to say about quantum chaos or L-funtions. 
The above discussion only serves as a background and a motivation for dicussing the 
probabilistic models that we will introduce. For a discussion of bounded variance 
in another context see PQ. 

2. Models and results 

2.1. The model. A point process on B C R with correlation functions p n (x%, . . . ,x n ), 
n > 1, has determinantal correlation functions if there is a function K : 

B x B — > R, the correlation kernel, such that 



(2.1) 



p n {xi, ...,x n ) = det(K(xi,Xj))? j=1 , 
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(2.3) u N {x)d N n{x) = ^ det{cb t (x 3 ))^ =1 det(^(xj))^id N p(x) 



n > 1. The interpretation of p n is that p n (xi, . . . , x n )dxi . . . dx n is the probabil- 
ity of finding particles in infinitesimal intervals dx\ 1 . . . , dx n around x\, . . .x n . In 
particular p\{x) = K[x, x) is the local density at x. 

Below we will construct kernels K by taking appropriate limits of other kernels 
and it is natural to ask if there is a determinantal point process whose correlation 
kernel is K. This can be answered using the following theorem. 

Theorem 2.1. Let Pk,N{%i, ■ ■ ■ ,%k)i k > 1, be the correlation functions of a de- 
terminantal point process on I C R with continuous correlation kernel K]y(x,y), 
N > 1. Assume that K^(x,y) — > K(x,y) uniformly on compact subsets of I 2 . 
Then there is a point process on I with correlation functions 

(2.2) p k (x u ...,x k ) = det(K( Xl .x 3 ))l j=1 , 

k>l. 

The theorem will be proved at the end of sect. 3. 

Let <f>i(t), tpi(t), i > 1 be functions in L 2 (X,p). Then, f\. [2lj. 

1 

~Z 

defines a measure on X N with determinantal correlation functions. If we have 
un{x) > for all x € X N and 

(2.4) Z=f dct(Mxj))%=ifet(Mxj))i!i=i dN rt x ) >0 ' 

Jx N 

we get a probability measure. We can think of a symmetric probability measure on 
X N as a point process on X with exactly N particles. In this paper we will have 
X = M or [0, oo) and p will be Lebesgue measure. The correlation kernel is given 

by 

N 

(2.5) K N (x,y) = J2 Mx)(A-%h(y), 

where A = (J x 4> l (x)%p j (x)dp(x))f, j=l . Note that det A = Z > 0. 

One natural way to obtain probability measures of the form l|2.3(l is from non- 
intersecting paths using the Karlin-McGregor theorem, |14|. Consider N one- 
dimensional Brownian motions started at y\ < ■ ■ ■ < ijn at time and conditioned 
to stop at z\ < ■ ■ ■ < zn at time S + T and not to intersect (coincide) in the whole 
time interval [0, S + T]. The induced measure on the positions x\ 1 . . . , xn at time 
S is then 

(2.6) Pn,s,t( x ) = \ det (Ps{Vi, Xj))" j= i det(p T (xi, Zj))? j=1 
where 

ft (x,y) = -i=e-^) 2 / 2t 
V 2nt 

is the transition kernel for one-dimensional Brownian motion. From the results 
discussed above it follows that (|2.6|) defines a point process on R with determinantal 
correlation functions. The correlation kernel is given by 

N 

(2-7) K n ,s,t( u > v ) = ^2 Ps(yk,u)(A~ 1 ) kj p T (v,z j ), 
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where 

(2.8) A = det(p s+T (y i ,z j ))^ 1 . 

In the limit T — > oo this model converges to Dyson's Brownian motion model, 
with (3 = 2, and was considered in see also [TT] . 

We can also consider the same type of measure but on [0, oo) and with appro- 
priate boundary conditions at the origin. We will consider reflecting or absorbing 
boundary conditions, where we have the transition kernels, [7], 

(2-9) P?(x,y) =Pt(x, y) +Pt(x,-y) 

and 

(2.10) pf(x,y)=p t (x,y)-p t (x,-y) 

respectively. We simply replace ps, Pt in 1|2.6|) with p™, p 1 ^ or p| b , p?^. In these 
cases we have initial points < y\ < • • • < y^ and final points < Z\ < ■ ■ ■ < z/y. 
We will be interested in these models as N — ► oo with fixed S, T or with S fixed 
and T — > oo. 



2.2. Correlation kernels. We want to obtain useful expressions for the correla- 
tion kernel 1|2.7[) with equidistant final positions, compare |13| . proposition 2.3. 

Theorem 2.2. Let Ti, : R 9 t — > L -\-it, L G R, and lei "/ be a simple closed curve 
that surrounds y\ < ■ ■ ■ < yi n -\; L is so large that Tl and 7 do not intersect. Set 
Zj = a(j — n) for some a > 0, 1 < j < 2n— 1. Consider the model \2. b)) with (yi) 
as initial conditions, and {zi)f" as final points, and with no boundary. Then 

ae ~(u 2 +v 2 )/2T f f 1 

(2.11) K2n+i,s,T{u,v) = —— T ———— / dw dz 







{2m) 2 S{S + T) J Tl 7 7 <*{w-z)t{T+S) _ 1 



2n-l p au,/(T+S) _ p ay 3 /{T+S) 

n 



e az/(T+S) _ e a Vj /(T+S) ' 

The theorem will be proved in sect. 3. 

The T — * 00 limit of this formula appears in J3|> compare theorem 12.31 below. 
There are analogues of the formula p.llfl for the absorbing and reflecting cases. 
We have not been able to write down a useful formula for the case of general final 
positions. The expression in (|2.11(l is more useful computationally than l|2.5|l but 
still rather complicated. We will obtain simpler formulas in certain special cases. 
First we will give a double contour integral formula for the case T = 00 in the 
absorbing and reflecting cases. We will also consider the N — > 00 formula in the 
absorbing case. Then we will specialize to the case when the initial points are also 
equidistant and T = 00 or T = S. In these last two cases we can obtain very nice 
formulas that are not in terms of contour integrals. The next theorem gives the 
analogue of proposition 2.3 in ^|| m the absorbing and reflecting cases. 



Theorem 2.3. Let Tl be as in theorem \2. 6 A and, assume that 7 surrounds y\, . . . , y^ 
and does not intersect Tl- Set Zj = j — 1 and assume < y\ < ■ • ■ < yjy. Then, 
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uniformly for (u,v) in a compact set in [0, oo) 2 , 
(2.12) 

lim K$ StT (u,v) = K$ s (u,v) 

1 — *oo 

N 2 2 



. 2C [ dw (dze^-v ) 2 /^ (e -(-^_ e -(^)>)J!!LfT , , 

(2tti) 2 5 y rj . y 7 w 2 - z 2 i = i z 2 _ y 2 



and 
(2.13) 

lim ^ SiT (u,w) = K% s (u,v) 



T^oo 

AT 2 2 

= / du; / dze(-« f {e -( z -uf /2S + e -(*+uf /2 S) ^_ TT 

(2tti) 2 S' y rt 7 7 w 2 - z 2 11 z 2 _ y 2 

The theorem will be proved in sect. 3. 

We will also write down a contour integral formula for the N — > oo limit of 
if^ s (u, v) under an assumption on the t/j-:s. (We could write a similar formula in 
the reflecting case, but we will omit it.) 

Theorem 2.4. Let Y l be as above, L arbitrary, and jm the two lines M 3 t — > 
=Ft ± iM wif/i M > 0. Le£ < j/i < j/2 < ■ ■ ■ assume that Yl'jLi ^/Vj < °°- 
Define 



(2.14) F(z) = I] \l , 

w/wc/i converges uniformly on all compact subsets of C. Set 

-, ,-L+Mi 

(2.16) K* S2 (u,v) = ^- e (w-vf/2S-( W -uf/2S dWj 

2™ J L — Mi 

and Kg(u,v) = K* si (u,v) + K* S2 (u, v). Then uniformly on compact subsets of 
[0,oo) 2 , 

(2.17) lim K$ s (u, v) = K a s b (u,v) = K* (u,v) - if* (-„,«) 
The theorem will be proved in section 3. 

We now come to the case when the initial points yj are equidistant. The next 
theorem is what makes it possible to compute the number variance in this case. To 
compute the number variance using the double contour integrals seems difficult. 

Theorem 2.5. Let yj = A + a(j -n), 1 < j < 2n - 1, < A < a, a > 0. Set 

(2.18) d= 2 -^. 

T/ien, uniformly for (u,v) in a compact subset ofM. 2 , 

(2.19) lim K N s {u, v) = K s {u - A,v - A), 

AT— >oo 
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where Kg(u,v) = a Lg(a u,a v), and 

1 , , p m{y+{2n-l)x) 

(2.20, Ls(x,v) = -ReZe-'<««-» nJ+ . (s _ xy 

Furthermore, if yj = aj, j > 1, then uniformly on compact subsets of [0, oo 
(2.21) lim K$ s (u, v) = Kf(u, v) = K s (u, v) - K s (-u, v), 



,2 



d 



an 

(2.22) lim K^ s (u,v) = K s e (u,v) = K s (u,v)+K s (-u,v). 

The theorem will be proved in sect. 3. 

The leading contribution to (|2.20|l comes from the terms n = and n = 1. The 
other terms are exponentially small in d. The leading part is 

sinTr(x-y) dcos7r(a; + y) + (y - x) sin7r(y + x) 

2.23) L s ^ ppT (x, y) = — + w , 

Tt{x — y) 7r(a z + (y — x) z ) 

so we have the ordinary sine kernel plus a non-translation invariant term. In par- 
ticular for the density we have 

^ 1 / cos(27ra/a) 

2.24 K s , appI (u, u = - 1 + V - 1 ' 

a \ ird 

so we have an oscillating density reflecting the initial configuration. 

It is also possible to give a more explicit formula in the case T = S. We have 
the following theorem. 



Theorem 2.6. Let 6i(x; u>) , i = 1,2,3,4, be the Jacobi theta functions, see 
let yj = A + a(j — n), 1 < j < In — 1, < A < a, a > 0, and let d be as in fJQgp . 
Then, uniformly on compact subsets of R, 

(2.25) lim K NtSjS (u,v)=K s ,s(u-A,v-A), 

N^oc 

where Kg t g(u,v) — a~~ 1 Lg t g(a~ 1 u,a~ 1 v) and 
(2.26) 

Ls,s{u,v) 



x 



2 (0;id)^6 3 {0:,id)9 4 (0;id) 

< 3 (u + v; 2id)9 1 (u - v; 2id) 9 2 (u + v; 2id)9 4 (u - v; 2id) 



sm 



h(7r(u - v) /2d) cosh(7r(w - v) /2d) 



If we neglect contributions which are exponentially small in d the leading part 
of ll2~2Tll is 

sinTr(u-u) cosn(u + v) 

(2.27) i5,S,appr(w, V) ~ 



2d sinh(7r(u - v) /2d) 2d cosh(7r(u - v) /2d) ' 

Note that as d — > oo both Lg and Lgg converge to the sine kernel. We see from 
(I2.23|l and i|2.27[l that Lgg decays much faster than Lg at long distances. 

We can also consider an averaged model by averaging over A in theorems 12.51 
and 12. 61 The averaged model has correlation functions 



(2.28) 



- / det{K( Xl - A,x 3 - A))? j=1 dA, 
a Jo 
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where K is the appropriate kernel Ks or Ks,s- This averaging will restore trans- 
lation invariance. The density in the averaged process will be a constant equal to 
1/a. In particular we have 

1 f a 

(2.29) - / K s (x- A,y- A)K s {y - A,x- A)dA 

a Jo 

sin 2 li^vl d 2_t^y\ 2 

a |_ V a / 

~ w 2 (x - y) 2 2w 2 a 2 (d 2 + {^f) 2 ' 
plus terms exponentially small in d. If we instead consider Ks,s we get 



(2.30) 



• 2 ir(x-y) -, 
sin — — 1 



Aa 2 d 2 sinh 2 8a 2 d 2 cosh 2 

lad lad 



2.3. The number variance. Let / C M be an interval and denote by the 
number of particles contained in /. We are interested in the variance, Var/f 
of this random variable in the determinantal point process with kernel K. The 
kernel K = Kn s,t that we have considered above is a reproducing kernel, i.e. 

(2.31) / K(x,y)K(y,z)dy = K(x,z). 



This is immeidiately clear from (|2.7|) and (|2.8|l and the same also holds in the 
absorbing and reflecting cases with integration over [0, oo) instead, and the repro- 
ducing property is inherited by the limiting kernels obtained above. Using (|2.31() 
and the determinantal form of the correlation functions it follows that 

(2.32) Var K {#I) = J dx J dyK{x, y)K(y, x) 

The sine kernel sm ^f~^^ a with density 1/a has the number variance, / = [0, L], 
(2.33) 

^sinckernelf-^) i> 
7T Z 



, 2ttL , 2wL.Tr n . ,2wL.. 2wL ,2wL s 

log + 7 + 1 + (- - Si ( )) - cos Ci ( , 

a a 2 a a a 



In the averaged models we get (we denote the averaging over A by < >), 

(2.34) < Vai K (#/)>= [ dx [ dy-f K (x - A, y - A)K (y - A, x - A)dA, 

J I J I- a J 

where K is Ks or Ks.s- The formulas above for the correlation functions and the 
formulas for the number variance are used to prove the next theorems. We will only 
consider the contributions from the leading parts of the kernels, H2.23|l and Ij2.27|l . 
Also, we will not use the reflecting and absorbing kernels. If the intervals are high 
up, I = [R,R + L] with R large, then the contribution from Ks(—u,v) in (|2.21|) 
and fl2~2l>|l will be small (like 0(1/ R)). 

Theorem 2.7. Consider the kernel Ks of theorem \2.5\ Define 9 and <fi by R/a — 
[R/a] = 8/w and L/a— [L/a] — <f>/w. The contribution to VarK s (#[R, R + L]) 
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coming from the leading part of the kernel I2.2tyl is, A — irL/a, 
(2.35) 

1 / cos 2(61 + (/>)+ cos 29\ (. 2nAd 



ird J \ y/A 2 + ir 2 d 2 

+ \ (l + 2A (| - Si(2A)j - cos 2A 

+ 7TT-, {(cos2(6» + (f>) + cos26)(h 3 (2iird) - h 3 (2A + 2md)) 
+ (sin2(6> + (/))- sin29){h 1 (2A + 2ind) + n- 2Si (2A))} 

+ IT-lr {M 2 ^ + 2i?rd) (sin 4(0 + <j))- sin46») - h 4 (2A + 2i7rd)(cos4(6> + <j>) + cos 49)} , 

07T U 

where h\(z) = 2Ref(z), h,2(z) — 2Imf(z), h 3 (z) = 2Reg(z), h 4 (z) = 2Img(z) and 



f°° sint , x f°° cost 

(2.36) f(z) = / — dt , g(z) = / —dt, 

Jq l -f- Z J Q I + Z 

for Rez > or Rez — 0, imz 7^ 0. These functions have the asymptotics f(z) — 
1/z - 2/z 3 + 0(1/ z A ), g{z) = l/z 2 + 0{l/z A ) as z -> 00. 

The theorem will be proved in sect. 4. 

If we average the expression for the variance over 9 or equivalently consider the 
averaged model, we get that the contribution to the number variance from the 
leading part is 
(2.37) 
1 



r2 



2-KLd/a „ ,2ttL, 2wL .it „ ,2ttL, 2ttL, 

log /r2 . = =+7 + l-Ci ) + ? -Si -cos ; 

^/L 2 /a 2 + d 2 a a 2 a a 



Apart from the logarithmic term we have exactly the same formula as for the sine 
kernel, l(2~35|> . It is not difficult to obtain (|2~37j> using the averaged correlation 
functions, i|2.29[) . The proof is then analogous to that in the sine kernel case, see 
sect. 4. The proof of (|2.34|) is a lengthy but rather straightforward computation. 

Note that when L is small compared to d the leading term in l|2.37|l is log(2%L/ a) 
which is what we have for the sine kernel. When L — > 00 the expression (|2.37|l con- 
verges to 

(2.38) ^(log(2^) + 7 +l) 

so the number variance saturates. Note that the saturation level does not depend 
directly on the mean spacing. If we rescale the model, see below, we get the same 
saturation level. 

We can also compute the number variance for the (S,S)-model where we have 
the kernel Ks,s- In this case we will only consider the averaged model. 



Theorem 2.8. Denote by Vd(L) the leading part of the averaged number variance 
of an interval of length L for the model with kernel Ks,s, i-e. \2.27jj , where d is 
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given by 12. Iffy . The 
(2.39) 



2 r L/a ( u/2ird \ 2 sin 2 w 2L f°° ( u/2ird \ 2 sin 2 



V d( L )- 7r 2j Q { S i n h(u/2ird)j u dU + na J„ L/a {sinh(u/2ird) > - 2<iu 
1 f°° min(u, L/2ad) 



u 



1 " cosh 2 u 



-du. 



Also, 



(2.40) lim — lim V d (L) = 1. 

d^oo -Aj log Ct L^oo 

If we compare l|2.39l) with the integrals which lead to Ij2.33|l , we see that there is 
a truncation effect which depends on d and which is responsible for the saturation. 
The limit 12.4U|) shows that the saturation level is similar to (I2.38|) for large d. 

As mentioned in the introduction the unitary group U(n) has been used in 
to model the (- function at height T, where n — \og(T/2ir). This n is obtained by 
equating the mean spacing in U(n), 2?t/n, with the mean spacing of the zeros at 
height T, which is 27r(log(T/27r)) _1 . Note that the eigenvalues of a random matrix 
from U (n) with respect to the Haar measure also show a kind of number variance 
saturation. The variance for the number of eigenvalues in an interval on T of length 
a, < a < 2ir is given by 

2 n n ~ 1 I 7 

, „ na na I v ^ n — k , ka 

(2.41) Var u{n) (a) = E^ sm Y 

k=l 

This increases as a function of a for < a < 7r and then decreases symmetrically, 
we have a maximum variance when a = w, i.e when we have a half circle. This 
maximum variance is 

(2.42) J_(i og (2n)+ 7 +l)+0(-), 

which is analogous to (|2.38|l if we set xd — n. Note that the number variance for the 
Riemann zeros at height T saturates at the mean level log(log(T/27r)) + const, 
[2], so equating the saturation levels (disregarding constant terms) leads to n = 
log(T/27r) again. This may be a more natural argument in a sense since it is not 
changed under rescaling. 

2.4. Approximation. Consider the situation in theorem |2.4l the absorbing (S, oo)- 
model, where the t/j's are not equally spaced but are given by yj — F^ 1 ^) for some 
increasing function F. If F is nice and do not vary too quickly, the t/j's will be 
almost equally spaced for long stretches of j, and hence we expect that the kernel 
(I2.17|l should be well approximated by l|2.21|) in a region where the average spacing 
is a. We will not attempt to make this clear in the greatest possible generality. Our 
goal is an approximation theorem valid for a certain class of functions F. Denote 
the correlation kernel Kg h (u, v) with initial points y = (j/j )5=i ^ ^s h (y' u ' v ) ^° 
indicate the dependence on y. We will prove the following approximation theorem: 

Theorem 2.9. Assume that F : [0, oo) — > [0, oo) is a C 2 -function that satisfies 

(i) F(x) < Cx 1+S , for all x > 0, for some 6 € (0, 1) and some constant C , 

(ii) F'(0) = and F'(x) > for x > 0, 
(Hi) F" is decreasing. 
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Fix a (large) and e > so that 1 + S + e < 2. Also, fix S e [F'(F(a), 1]. Define 
V = {Vj)%i °y Vj = F^U), 3 > 1 V-j = ~Vj f° r 3 > 0- For m > 1 set 

ir, in\ ^ _ tym — Vm—j — Vm+j 

1 ' Km ~ ^Tu -v -)(v ^ -v )' 

j—\ \y m ym—j J\ym+j ym) 

Ki 1 = F'iVm) and rj m = F"(y m ). Also let £ m = y m - £ m S. There is an m = 
m(a) > 1 such that |£ m — a\ < X m . We have roughly m(a) ~ F(a). Set A(a) — 
A m(Q) ~ (F'iF-^a))- 1 , 77(a) = n m{a) ~ ^"(^(o-)) and 

(2.44) T (a) - min( — i=, mCcO^'Wi-*)). 

Define y = 67/ j/j = A(a)j ; j > 1. vlsswme i/iai < T < To (a). There are 

constants Co, C , which depend on F and e, but not a, such that 

(2.45) | e -«°0(«-« )tf u , „) _ ^ b (y; u - C(«), « - C(«)) 

< C[A(a)S- 3 / 2 e- fl2 /ss + (T 2 + ^^-(i-^/a+^-i^ 

/or all u,v & [a — T, a + T] provided R 2 lies in the interval 

O \ 2 /(!- e ) Til + £Q N 



(2.46) 



c n max 



,A(a 



Note that Kg h (y;x,y) is given by (|2.21|l with a = A(a), so provided the right 
hand side of i|2.45[l is small we have an approximation with a kernel having equally 
spaced initial points. The factor exp(— £(a)(u — v)) in front of Kg h (y; u, v) in (|2.45(1 
does not affect the correlation functions corresponding to this kernel since it cancels 
in the determinant. Let us consider two examples of theorem 2.8. 

Example 2.10. Let F(x) = x 1+s with < S < 1 and fix e > small. Then 

X(a)- 1 ~ (1 + S)a s ^ 1+S \ 77(a) ~ 5(1 + S)a-^- s ^<- 1+ ^ and m(a) ~ a 1+s as 

a -> 00. Choose 5 = 1. If < T < a (1 -^ 2 / 2 , then 

(2.47) 

K ab (| ,. u) _ ^ab (i; u _ c(a)) „ _ c(a)) | < 

for all u, v € [a — T, a + T]. We see that this is only interesting if 6 < 2 — \/3, 
otherwise the right hand side of (|2.47|) does not go to zero as a — ► 00 unless we also 
let T — > 0. It also follows that 

(2.48) km ^MW^a+^a + Xj^liifH 

a^oo - A(a) A(a) 7r(x — y) 

uniformly for x,y in a compact set, so as we go up the line we see the sine kernel 
process. Note that when x and y belong to a compact set we can take T = CX(a) 
and the right hand side of (|2.47(l goes to zero as a — ► 00 for S < 1 provided we choose 
e small enough. Thus we can extend l|2.48|l to all < 5 < 1 We see a transition 
from a non-universal regime for small a to the universal sine kernel regime for large 
a. 

Example 2.11. Let F(x) = ^ log #- — #-. This F does not satisfy all the con- 
ditions in theorem 12.91 but we can modify it for small x so that it does without 
changing the yy.s. Then A(a)j^log^- — > 1, 2irarj(a) — ► 1 and m(a)/F(a) —> 1 as 
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ot — ► oo. Fix eo > small and let 5 = e where (1 — e) 2 /(l + e) — 2e = 1 — eo. Pick 
R = coT^ 2+e a e . The theorem then shows that if 1 < T < a^ 2 ~ e , then (S = 1) 



(2.49) e -««-«) K x (y; u, v) - K^y; u, v) 



< CT 2 a~ 1+e 



for u 1 v € [a — T 1 a + T] . We see that the parameter d for the aproximate kernel is 

2tt 8tt 3 



d(a) = 



A(a) 2 (log(a/27r)) 2 ' 
The saturation level is the (disregarding constants) 

(2.50) ~ \ log d(a) ~ A lo g( lo S T")- 

This differs by a factor 2 from what we would like to have if we want to model 
the ^-function. Note that the mean spacing at height a is the same as for the £- 
function. However the saturation level does not change in the equidistant case if we 
rescale the point process linearly (the distance a between the yj'.s is not changed, it 
is the constructed point process that is rescaled) . Hence the mean spacing and the 
saturation level are independent. We can make a better model of the ^-function by 
picking a suitable F, see below, and then make a non-linear rescaling. Below we 
will construct a model for the unfolded zeros Sj = Af(Ej). 

Consider a determinantal point process on [0,oo) with correlation kernel K and 
let G(x) be a strictly increasing C 1 -function. Then 



(2.51) K(x, y) = K(G(x), G(y)W G> (x)G> (y) 

defines a new, rescaled point process. The density for the new process is p(G(x))G' (x) 
if p{x) is the density for the original process. 

Let yj — F^ 1 ^), j > 1, where F satisfies the conditions of theorem l2.9l The den- 
sity at x is then F'(x) and if we rescale with G we get F' (G(x))G' (x) = -^F(G{x)). 
Hence to get a constant density we should pick G(x) = F^ 1 (x). Choose 

(2.52) F{x) = 1+ T { i og ±)U2 dt 

for x > 27re and define it for < x < 2ire so that the conditions in theorem 12. 91 are 
satisfied. Consider the model corresponding to this F and rescale it using G = F^ 1 
as in (|2.51|) . We will say "approximately" below without beeing too precise. The 
estimates involved can be made precise with a little effort. By theorem !2.9l at height 
T, the kernel will be approximately (5 = 1) 



(2.53) K^F-^x) - F-\T), F'^y) - F ^(T)) ^(F^y(x)(F-^'(y) 

for x,y in a neighbourhood of T, where yj = X(F~ 1 (T))j , j > 1, A(a) _1 ~ F'(a). 
(We have approximated C,(a) with a.) Now, for x close to T we have F^ 1 (x) — 
F- X (T) w (x - TXF-^'lT). Set b = (F _1 )'(T) w A(T). By ll2~53l and $T%% at 
height T (large) the kernel will approximately equal 
(2.54) 

T ( _T - sin7r(a; - y) | rfcos7r(x + y - 2T) + (y - x) sin7r(a; + y - 2T) 

■k{x — y) -K{d z + (y — xY) 

where d = 2tt/A(T) 2 . We have A(T) - (log^)" 1 / 2 . As T -> oo , d oo and 
we see that the kernel in l|2.54ll converges to the sine kernel. The point process we 
have constructed thus has the correct universal asymptotics as we go up the line, 



12 



K. JOHANSSON 



and it is non-universal for small T, since it depends on the particular F we have 
chosen. The saturation level for the kernel in (|2.54(l is 



(2.55) — log(27rd)~ — log(log— ) 



1 . ,„ „ 1 . T 
2n' 

when T is large, which is exactly what we would like to have at height T. We 
have thus constructed a determinantal point process in [0, oo) which has, in some 
aspects, similar behaviour to the unfolded zeros of the ^-function. 

2.5. Correlation kernels close to the origin. Consider the kernels K$ on R, 
and Kg h and K™ on [0, oo) in theorem l2.5l If we are interested in say the distribu- 
tion of the first particle to the right of the origin we have to compute the probability 
of having no particle in [0, £] . Let X be the position of the first particle to the right 
of the origin. Then, if the correlation kernel is K, 

(2.56) ¥[X < f] = 1 - P[no particle in [0, £]] 

= 1 - det(7 - K) L 2([ ,Q), 

where the second equality is a standard result for determinantal point processes. If 
we vary a and S in such a way that S/a 2 — > oo, it follows from theorem 12 . 51 that 

, n t ts i \ sin7r(u — v) 

(2.57) lim aKslau, av) = : r — , 

5/a 2 ^oo TT(U — V) 

a K/ ^ sin7r(u — v) sin7r(u + u) 

(2.58) lim aKf{au,av) = , \ 

S/a 2 ^oc TT(U — V) TT{U + V) 

and 

,„ ,„> ,. j-^ro/ \ sin7r(M-w) sinir(u + v) 

(2.59) hm aK>s c (au, av) = ^ + v / 

S/a 2 ^oo 1T{U — V) TT[U + V) 

These kernels can also be obtained from the classical compact groups and have been 
used by Katz and Sarnak to model the lowest zeros in families of L-functions, see 
|15j . The above results show that they can also be otained in a natural way from 
non-intersecting paths with appropriate boundary conditions. 

The kernels in the right hand side of Ij2.58|l and Ij2.59|l are directly related to 
special instances of the Bessel kernel, 

(2.60) *,(„, „) . ''"MW.W-'fWJ^W, 

where J u (x) is the ordinary Bessel function. In fact a simple computation shows 
that if we define the rescaled Bessel kernel B v by 

(2.61) B v (x, y) = V2TT 2 x2TT 2 yB l/ {n 2 x 2 ,n 2 y 2 ) 7 
then 

„„. ~ . , sin7r(M — v) shnr(u + v) 

2.62 B ±1/2 (x,y ) = —A r 1 T ( \ / • 

' ir(u — v) tt{u + v) 

When v is an integer the kernel B v appears in the scaling limit around the smallest 
eigenvalue in LUE, the Laguerre Unitary Ensemble. If x\, . . . , xn are the eigen- 
values of M*M, where M is a (y + N) x N complex matrix with independent 
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standard complex Gaussian elements, N(Q, 1/2) + iN(0, 1/2), then xi, . . . ,xn is a 
finite determinantal point process with correlation kernel K 1 ^ and 

This interpretation does not work for v = ±1/2. 

3. Computation of the correlation functions 

In this section we will use the formula (|2.7|) to compute the correlation functions. 
If A is a matrix and b a column vector we will denote by (A|6)fc the matrix where 
column k in A is replaced by b. Let 

p= (pt(v, Zq) . . .pr(v, z N )) T ', 

where Zj = a(j — n), < j < N — In. By 1|2.7|) and Kramers rule we have 

det(A|p) fc 



(3-1) K 2n ,s,T(y;u,v) = ^2ps(yk,- 



k=0 



detA 



where A = {ps+rilJi, Zk))j n k=0 - If C M is the contour f t + iM, M G K, t G K, we 
have 

(3-2) Pr(^^ = \f^^ v2/2T J c e^ 2 /Wfe(r),^)e Sfc(T)2/2(T+5) dr, 
where y fe = jf fe (r) = i>(T + 5)/T + iry/S(T + S)/T; set also y.^y.iii^k. Then, 
(3-3) 

det(A|p) fc _ j T + S ^ v 2 /2T r c _ T 2 /2+ijk(T) 2 /2(T+s) det(pT+s(yj, Zi))% =0 
detA V 2vrT 7 Cm det(p T+s ( % , 

Since = a(j — n) the determinants in the quotient in the right hand side of <|3.3|) 
can be computed using Vandcrmonde's determinant and we find 



„ay k /(T+S) „a%n/(T+S) 
e ay k /(T+S) _ e a yj /(T+S) ' 



,(wJ-{/2)/2(T+5)+on(y*-5 fc )/2(T+S) TT 



Inserting this into l|3.2[) and making the change of variables w = jjk (t) we obtain 

/g 4 x det(A|p)fc _ 1 c -t, 2 /2T+^/2(T+S)+» ? , t a/(T+S) 

detA iV^u? 



/ d u;e T(«'-(T+S) 1 ,/T) 2 /2S(T+S)-n»a/(T+S) TT £_ 



aw/{T+S) _ p a Vj /{T+S) 



Vk/(T+S) _ e a V] /(T+S)- 



We can now use the expression in l|3.4|) and insert it into (3.1) to get 
(3.5) 

2n 

K 2niS ,T(y,U,v) = J_ e -v 2 /2T+u*/2TY e -T(y k -(T+S)u/Tf/2S{T+S) 
— v.tt ,S ^ — ^ 



2nS 

k=0 



f a 

x / dwe T{w ~ {T+s)v/T]2 /2S{T+S) e na{vk ~ w)/{T+s) TT — 
Jr, ?A e a 



aw/(T+S) _ e a Vj /(T+S) 



Vk/{T+S) _ e a Vj /(T+S) ■ 
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This is the basic formula from which the others will be derived. It is now straight- 
forward to prove theorem 12. 21 

Proof, of theorem \2.2\ If we apply the residue theorem in l|2.11|l we get the ex- 
pression l|3.5[) . Since and 7 do not intersect the w = z singularity does not 
contribute. □ 

We turn next to theorem l2.3l 

Proof, of theorem \2.!ft Consider first the absorbing case. We use the formulas (|3.1|) 
and (|3.3|1 but with pf 3 instead of pt- The evaluation of the two determinants can 
now be done using the following Vandermonde type identity 
(3.6) 

N 

det(e (2 ^ 1) ^-e- (2l - 1) ^)i V J =i =n( e3; '- e ""') U ((^-e^j'-^'-e^f)- 

i=l l<i<j<N 

Using this identity we find 
(3.7) 



det( A\p ah ) k IT + S y't f ^e 

— ■ - e 2(T + S) 2T I e~~ 



T+*\/ T(T + S) „ T i\J T(T + S) 



detA V 2nT ~ J Cm ~ - e -ft- 



x 



n 



e T+ l \j T(T + S) _ g T l \l T(T+S) j _ /g^fe _ ^ 



(e T 



jjtk e T + s — e T + s — e T + s — e T + s 



In this expression we can expand the exponentials and take the T — > oo limit. We 
find that 

(3.8) hm d6t ^> = ' [ e W«"» TT 

uniformly for (u, v) in a compact set. Here we have changed integration variable 
by putting w = v + ir\/~S. Hence 

(3.9) 

lim K^ T (u,v) 

1 — >oo 

N 

> . ,\ / • . . </• ^— - ir~ - ir 

-dw. 



1 V ( -{y h -uf/2S -(y k +u) 2 /2s\ [ Jw-vf/2S^_ TT 



,2 „,2 



k=l Vk j7 tk y k ~ y j 

That this expression equals the expression in (|2.12() follows from the residue theo- 
rem. 

The proof of (|2.13|) is completely analogous. Instead of (|3.6|) we use the identity 
(3.10) det (e^- 1 ^ - e-y-^H = J] ((e^ + e"^) - {e Xi - e~ Xi ). 

h0= l<i<j<N 

a 
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Proof, (of theorem \2.5} . Set 



A / 2 " 

2T 



(3.11) F N (z) = Y[ 1-- 



i=i\ ^' 



so that, by (EHJ, 
(3.12) 



1 f , f > 2 / (*- ;! ) 2 i£±^\ 2w Fjv(u>) 

I I J ~ „ <T3 I * ^ ^ 1 



V) = (2^5 j To dw 4 r " £ , - , 2 F„ W 

Here 7 M is the curve given by t — ► — i ± im, i > r and s — ► r is, — Af < s < M, 
where < r < y\. The fact that F/v(z) — > -F(z) uniformly on compact subsets of C, 
together with estimates like (|5.10|) and (|5.11() below, which can be used to restrict 
the z— and w— integrations, shows that 



(3.13) lim Jf»^( UjV ) = 7 __- / dw / dze^ 



x e^r~ - e~^~ — ^-prH- = K% b {u,v). 

\ J w 2 - z 2 F(z) 6 v ' ' 

uniformly for (u,v) in a compact set. Let 7^ r be the contour which is the image 
of 7 Mr under z — » — z. Using 

2w w ( 1 1 



Ur — z^ z \ u> — z w + z 

we see that Kg h (u, v) = K s (u, v) ~ K s (—u, v), where 

r „. , 1 [ , [ , > 2 ('-" > 2 1 tuFfio) 
K sM = -— T - dw dze^z ^ 



(2ttz) 2 5' 7 ro 7 7M r " z-w zF(z) ' 

7m, r = 7m r + Tajt- ^ we ^ r ~ * 0+ we pick up a contribution from the pole 
z = w. This leads to K s = Kg 1 + i"Tg 2 as stated in the theorem. □ 

Next we consider theorem 12. 51 

Proof, of iheorem \2.5\ It suffices to consider the case A = 0, otherwise we replace 
(u, v) with (u — A, v — A). From the proof of proposition 2.3 in |13| . the T — > 00 
limit of (|3.5|l is 

(3.14) K NtS (u,v)= lim K n , s ,t(u,v) 

T — >oo 

1 x (ak-u) 2 / ( m -f) 2 -r-r w — a j 



Now, 



27ri5' /r, - LJ - ak — aj 

k=-n LL j=-n,j^k J 



/ -1 \ h lb / > \ ^ — 1 

w — aj [—1)W "I — r / 10 \ -p-r n — J 



(3.15) TT 4Zf2 7 = iZi^TT 1- ™ n , 
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It follows from (|3.14|) and (|3.15|l that 
(3.16) lim K NS (u,v) = K s (u,v) 

N— >oo 

uniformly for (u, v) in a compact set. To prove the convergence we need some 
estimates so that we can cut off the fc-summation and the w-integration. We omit 
the details. Now, 

OO / 9 \ 

3.17 w TT 1 -_— =_ sin -w. 

fj[\ W) J T a 

Hence, 

where L is arbitrary. Replace w by u> + v — u + ak in the integral in l|3.18fl and use 
Cauchy's theorem to get 



OO 

a 



(3.19) Ks (u,v) = e^^e^^±- '-dw 



2lT 2 iS T^q Jv -i W + V — U 

a f u, 2 ~2 UW kan, sin Z-(w + v — u) 
+ 7TT-^ > / e ^ e— ^ -dw 

l-K lb , z — ' Jr. W + V — U 



h— — oo 

e { W 2 -2uwD)/2S gin I (y, + w _ u ) 



dw, 



2ir 2 iS J r e aw / s - 1 id + u - u 

where V = T-\ +Y\. The function e aw / s — 1 has simple zeros at w = 2TrSni/a, 
n E Z and hence by the residue theorem applied to the last integral in 13.19fl . 

3.20 K s (u,v) = -Ve /Q - 27rtira / Q -■ 

If we set L(x, y) = aK(ax, ay) and define d as in (|2.18[1 we find 

(3.21) L(x, y ),iy e -irdn»-2 TOfa Bin ff (nd» + t/ -x) 

1 „ — 7rdn 2 — 7rnd 

_ _L e 7ri(i/-(2n+i)x) 

27r ^— ' —rz.d + i(u — a;) 

_ _ _ e 7ri(y+(2n-l) a; )^ 

2tt ^ —nd + i(y — x) 

If we change n to — n in the first sum we get (I2.2U|) . 

In the absorbing case we have by (|3.9() . (|3 - 1 1 1) and the identity 

w -pr w 2 - y] = f_l 1_\ F N {w) 

Vk j-^ yl-Vj \w-y k w + y k J F' N (y k ) ' 
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that 



OO r. 

(3.22) ^^K«)~E/ r 



1 1 



e 23 — e 53 



(— l) fc sin —wdw. 



w ~yk w + y k J a 

If we now note that 

(afc-t.) 2 (ak + u) 2 \ / 1 1 



X , , / (afc-u) 2 (afc + u) 2 \ 

(3.23) £(-l) fe e-^^ - e"^^ 

fc=i ^ ' 



io - yh w + y k 



(.*-.,)' (-l) fe V-^ (.* + .,)' (-l) fe 

= \ e 2^ — — I \ e - — — 

w — ak ^— ' w — ak 

fcez feez 

then (|3.18l) and l|3.22[) give (|2.21() . The proof of (|2.22l) is completely analogous. □ 

The proof of theorem l2.6l is a similar but somewhat more complicated computa- 
tion where we have to use 0-function identities. 

Proof, of theorem \2.6l Our starting point is the formula i|3.5|) with T = S. We can 
assume that A = so that yj = a(j — n), < j < 2n. Set 

n n 

(3.24) F n (z; q) = \{(\ - j z) J](l - <?/z), 
where 

(3.25) q = e - a2 ' 2S = e-* /d . 
Then (|3.5|l and a short computation gives 

ae _(. u 2 _ n 2 )/2 S n n 1 

(3.26) K 2n+1<s ,s(u, v) = 2{2m)2s2 dw ^ dz ea{w _ z)/2S _ i 

\(w-2v) 2 -(z-2u) 2 ]/2S F n (e aw/2S j g) 

F„(e-/ 2S ;g)' 

where 7 surrounds = a(j — n.), < j < 2n and does not interesect T^. Let 7m 
be t — > ^ft ± iM. Replacing 7 by m (|3.26() we pick up a contribution from the 
pole z = w, when |Imw| < M. Hence 

ae -( v 2 -u 2 )/2S , , t 

(3.27) K 2n +i,s,s(u,v) = dw dz 



2(2ni) 2 S 2 J Tl J 1m e a ( w ~ z )/ 2S -l 
F n (e a */ 2S ;q) 

e i(w~2v) 2 -(z~2u) 2 ]/4S dw _ 

L — iM 



2ni 



In this expression we can control the n — > 00 limit (using some estimates of F n , 
compare (I5.1U|) . I|5.11[) ). Set 

OO OO 

(3.28) F(z-q) = [](1 - tfz) - q j /z). 

3=0 j=i 
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It follows from (|3~2T)) that 
(3.29) 



2 -« 2 )/2S r r i 



Ks,s(u,v)= lim K 2 n+i,s,s{u,v) = ————— — I dw I dz 



at^oo "»+w>\ ' > 2{2mYS 2 J Tl J lM e a(«-»)/as _ i 

F(e«*/ M ; g) 

- / e [(«--2 t ,) 2 -( 2 -2«) 2 ]/4S du; ^ 

We now compute the z-integral in l|3.29[l using the residue theorem. Apart from 
the pole z — w if |Imu>| < M we have simple poles at z = ak, k G Z. We obtain 



[(iu-2u) 2 -(2-2m) 2 ]/4S 



e -(i, 2 -« 2 )/2S f 

(3.30) K s ,s(v,v) = — V/ dwe 

1 F{q~ w / a ) 
q-w/a+k q-kp/fq-ky 

In the integral in l|3.30|l we make the change of variables w — > a(w + fc) and use 
Cauchy's theorem. A computation shows that 



(_1\k-l 00 

(3-32) 9 - fe F'( g - fc )= y^nt 1 ^") 2 

n— 1 



which gives 



If we set 



1 ^""'-na-o-^n,-"). 



fl-™ - 1 q- k F'(q~ k ) ~ 1 ' 1-q 



(3.33) G(z;g) = JJ(l-^)(l-9»7«), 

i=i 

we obtain 

np -{v 2 -u 2 )/2S r 

(3.34) gg, g (u,t>) = n 5 V/ dwe A[(™+^ 2 «) 2 -(- 



k-2u) 2 } 



y , fc a 2 „/25 g(g" W ;g) 

G(i;<z) ' 

Make the change of variables u? — > uj + (t> — u)/a and perform the fc-summation to 
get 

ae -(u 2 +« 2 )/2S /■ e (aw-u-v) 2 /S Qfq-w+(u-v)/a. g\ 

(3.35) = — s I e _ Batu/s _ 1 G(1;g) ' dw, 
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where T = Ti — r_i as before. The integrand in (|3.35|1 has simple poles at w = dni, 
n € Z, and the residue theorem gives 

p -(v 2 +u 2 )/2S „ fi( n -dni+{u— v)/a. n \ 

(3.36) K St sM = <£ e ^»-vf/s G{q , <?) 

A computation leads to 

(3.37) L s ,s{u, v) = aK s , s (au, av) 

-Tr(u-v) 2 /2d\^ -Trdn 2 /2-7rni(u+v) G , ((-l)"e~ 7r ^~"^ d ; q) 

^ G(l;q) 

If we write q = e™ , uj = i/d, p = e™, C Q = Un^i 1 ~ <? 2 ™)' [t follows from the 
product representations of the ^-functions that 

(3.38) G(p*;q)- 1 Wxiufotev) 



G(-p 2 ;q) 



2C^q 1 / i sin the 
1 9 2 (x; u>)9 3 {x; u) 



2C$q 1 / 4 cos 7tx 

Write u>' = —2/u> = 2id. If we insert the formulas (|3.38|l into (|3.37|l we obtain 
(3.39) 

T<(V,-V) 2 

Uj 2CggV4 G (l; q) \ e sin - (u _ u) 

, c w(»-i) 2 -^(2n-i) i (u+.) g 2 (f ("-^);^)g 3 (f ("-■"); ^) 

cos^(u-v) 

- ■— ' %(u + v; <J)9 X {% (u - u); w)0 4 (f(u - «); w) 



2C V /4 G(l;g) sin 2^( M _ w ) 
2 (u + u; o/)0 2 (f (u - «);a;)03(S(u - «); w) ' 



cosf^u- 

where we have used the series expansions of the ^-functions, see (I3.43|) below. 

We can simplify l|3.39|l somewhat by using some 9- function identities. The first 
is Jacobi's transformations: 

a , \ 9 k (x;~l/uj) 
9 k {ojx;uj) - 



k = 1, 2, 3, 4, where a\ = —i and a 2 — a 3 = 0:4 = 1, and the formulas |17j . p. 17, 

(3.40) 9 1 (x;lu)9 2 (x;lu) = 9 1 (2x; 2w)6> 4 (0; 2oj), 

6 3 {x; uj)9 4 (x; uj) = 9 4 {2x; 2cj)6» 4 (0; 2u>). 

This leads to 
(3.41) 

i9 i {0;~2/uj) \ O^u + v;uj')9 1 {u-v-u') 9 2 {u + v; u')6 4 (u - v; uj') 
S{U,V) ~ 2wC 2 q 1 /4 G (i ;(? ) { _ zsin ^( u _ w ) + - icos 7f( u - v ) 
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Now, by the product formulas for the values of the 0-functions at the origin we have 
G(1|?, -ll (1 ' } _ 2nC%qV* 

and consequently 

(3.42) 2C V /4 G(l;g) = -^(0^)^(0;^) 

= -^#4(0; --)# 3 (0; --)# 2 (0; --) 2 

ui 1 uj to to 

by Jacobi's transformation and the formula 
By Landen's transformation 



4 (O; -2/w) = ^3(0;-lM(0;-lH 

and hence 

i0 4 (O;-2/w) _ 



2wC V /4 G(l;<z) 2 (O; -l/w)V^(0; -1/^(0; -1/w) ' 

If we insert this into (|3.27|) we obtain 12.2611 and the theorem is proved. That the 
leading behaviour of the kernel is given by 12.27fl follows from the series expansions 
of the 9- functions: 

(3.43) 



»i(sc; 


r) 


_ j ~y^^_^n e iriT(n~l/2) 2 +Tri(2n 






nGZ 




r) 


_ ^""^ e TTiT(7i-l/2) 2 +ni(2n-l)x 








k(x; 


r) 


\ ^ -nirn 2 +2irinx 






neZ 


k(x; 


r) 


^ V ^TTtrn 2 -\-27rinx 









□ 



We should also prove theorem 12. II 



Proof. We will use the following criterion of Lenard, [T^], |23| : 

A family of localy integrable functions pk ■ I k — > K, k = 1, 2, . . . are the correla- 
tion functions of some point process if and only if the following two conditions are 
satisfied 

a) (Symmetry) For any a S Sk, 

Pk{x a (i), ■ ■ ■ ,x a (k)) = Pk(xi, ■ ■ - ,x k ). 

b) (Positivity) For any finite set of measurable bounded functions <fik ■ I k — ► K, 
k = 0, . . . , N, with compact support, such that 

N 

(3.44) &> + X) <t>k{x ix ,...,x ih )>Q 

fe=l i 1 ^---^i k 
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for all (x%, . . . , xn) € I N it holds that 

N 

(3.45) ^o + J^/ 4>k(xi,...,Xk)pk(xi,...,Xk)dxi...dxk > 0. 

k=i Jl " 

The uniform convergence of to K on compact sets implies that K is continuous 
and hence pt is locally integrable. It is also symmetric. We know that l|3.45fl holds 
with pk t N instead of pk since pk,N are the correlation functions of a point process. 
Since all 4>k have compact support and are bounded we can take use the uniform 
convergence of Kjy to K and take N — > oo to get (13.45(1 . This completes the 
proof. □ 



4. Computation of the number variance 

We will first show how l|2.37f) can be obtained from l|2.29f) . By (|2.34|) we want 
to compute 

/>L />oo />L />0 

(4.1) dx dyf{x-y) + dx dyf(x-y), 

Jo Jl Jo J -oo 

where 



(42) f(x) Sln2 g I S ~ {x/a)2 

K ' J{ ' ir 2 x 2 2ir 2 a 2 (d 2 + {x/a) 2 ) 2 ' 

Since f(x) is even we see that H4.ll equals 

/0 />oo />L />oo 

dx / dyf(x -y) = 2 / dyf(x + y) 
-L Jo Jo Jo 

pL / poo \ pL poo 

= 2 / ( / ) dx = 2 / xf(x)dx + 2L f{x)dx. 

Jo \Jx J Jo Jl 

If we take f(x) = (sin 2 irx)/n 2 x 2 in 14. 3|) we get l|2.33|l using (|4.16|l and (14.171 
below. Inserting l|4.2fl into l|4.3fl and computing the integrals we obtain (|2.37|) . A 
similar computation using l|2.29|l leads to l|2.39fl . To prove (|2.4U|) we have to show 
that 

iAA . 2 f°° u 2 sin 2 27rud J 

(4.4) hm / ^ du = 1. 

d^oologdJo sinh u u 

Using sinh 2 u > u 2 we get 

2m 2 sin 2 27rud, f 1 1 - cos 2nud , f°° du 



(4.5) / - — 2 du < / du + 



sinh 2 u u J u J 1 sinh 2 

w 1-cosj;, f°° du 
dx 



x J i sinh 2 u 



= log(2vrd) + 7 - Ci (2?rd) 

and hence 



du 



sinh 2 u ' 



2 Z" 00 u 2 sin 2 2nud , 

limsup- ; / 5 <xu < 1. 

d^oo log d Jo sinh u u 
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Given e > we can choose S > so that \x 1 sinhrc — 1| < e if |x| < 6. Thus 

, x f 00 2u 2 sin 2 2TTud , 2 f s sin 2 2nud , 

(4.6) / 2 cfc > ^ / cZm 

Jo sinh u u [1 + J Q m 



l -^{\og{2TT d5) + 7 - Ci (2tt^)), 



which gives 



(l + «0 



,. . , 2 Z" 00 u 2 sin 2 27rwd , 

limini : / ^ au > 1, 

d^oo logd J sinh m u 



and we have proved 14.4JI . This completes the proof of theorem l2.8l 
Proof, of theorem \2. 7| We have 

rR+L roo 

Vav Ks {#[R,R + L})= dx dyK{x,y)K(y,x) 

JO JR+L 



rR+L pR. 

+ dx dyK(x,y)K(y,x). 

JR J-oo 

Set £*(*,!,) = iL(f,£), and 

v(A,9) = dx dyL*{~x + 9,y + 9)L*(y + 9,-x + ( 
Jo Jo 

A computation shows that 

(4.7) Vwc Ks (#[R, R + L}) = v{—,6 + 0) + v(—,-6). 

a a 

Hence we have to compute v(A, 9). If we set 

e -2i({n+k)x+(m+k)y) 



(4 - 8) G ^ k ^J Q dx J Q d y^ dn + i(y + xm,dm + % + ,))' 

then 
(4.9) 

u ( y l ) 6)) = _L e - dn{n -^e 2l{n - m)e e- dm{m+1) G(n,m,0)+e- dm{m - 1) G{n,m 7 ~l) 

If we neglect the terms that are exponentially small in d we get 
(4.10) 

An 2 v{A, 9) = -G(0, 0, 0) + G(0, 0, -1) + G(l, 1, -1) + e 2ie [G(l, 0, -1) - G(0, -1, 0) 

- G(l, 0, 0)] + e- 2l0 G(0, 1, -1) - e 4l9 G(l, -1, 0) + c.c. 

(Note that G(0,0, 0) and G(0,0, — 1) are not individually convergent but have to 
be considered together.) 

We will now outline how (|4.1l)|) can be computed without giving all details. Set 

rA ^ -2mi 

(4.11) H 1 {A-n,m)= — — dx, 

J irdm + ix 

A> 0, n,m € Z, d> and 

(4.12) H 2 (A;n,m) = / — -dx, 

J 4 irdm + tx 



-c.c. 
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A > 0, d > 0, n 7^ 0, 7i, m € Z. Some computation now gives 
(4.13) 

^ ^ — 2«(n— m)A 

G(n, m, k) = . — \H 2 (A\ m + k,m) — H 2 (A; to + k, nil 

Imdyn — m) z 

' [ifi (A; m + k, n) - H x (A; m + k, to) + H 1 (A;n + k,m) - H 1 (A;n + k,n) 



2md(n — to) 2 
if n =/= to, 
(4.14) 

G(n, n,k) = l- e - 2l(n+k ^ A - i log (A 2 + ir 2 d 2 ) + % arctan ^ + logfrd) - ijsgn (n) 

Z ./l Z 

- 2(4 - nidn)(n + k)H 2 (A;n + k,n)- 2mdn(n + k)H 2 (0;n + k,n) + iH 1 (A;n + k,n) 
if n =^ and 

(4.15) G(0, 0, jfe) - G(0, 0, 0) = 1 - e" 2iM - 2fcAff 2 (A; fc, 0) + H X (A; k, 0). 

The next step is to express the functions H x and H 2 in terms of f(z) and (j(z) 
defined by l|2.36[l . and in terms of the sine and cosine integrals: 

(4.16) Si(A)=f A ^y d y ; \ - Si (A) = f ^dy, 

Jo V 2 J A V 

pA -t poo 

(4.17) / ~ C08V dy = 7 + log A- Ci(A) : Ci(A) = - ^-dy. 
Jo V J A V 

After some computation we obtain 

(4.18) H X (A; n, to) = — i log \J A 2 + Tt 2 d 2 m 2 + i log |7rdm| + arg(4 — mdm) 
+ &rg(iridm) + iG(— 2iridm\n\) — i(cos 2nA)g(2A\n\ — 2iridm\n\) 

+ i(s'm 2\n\A)f(2A\n\ - 2iridm\n\) + f (-2iridm\n\)sgn (n) 

- (cos2nA)/(2,4|n| - 27ridTO,|n|)sgn (n) - (sin 2nA)g(2A\n\ - 2iridm\n\), 

if n =/= 0, to =/= 0, 

(4.19) n, 0) = -i(7 + log(2A|n|) - Ci (2\n\A) - sgn (n)Si (2A|n|), 
(4.20) 

H 2 (A;n,m) = - i (cos 2n A) g (2 A\n\ - 2iridm\n\) + i (sin 2\n\ A) f (2 A\n\ - 2iridm\n\) 

- (cos2nA)g/(2A|n| - 27r«dTO|n|)sgn (n) - (sin 2nA)g(2A\n\ -2mdm\n\), 

if to ^ 0, n ^ 0, 

(4.21) # 2 (A; n, 0) = iCi (2A|n|) - sgn - Si (2A|n|)) 

if 7^ 0, and finally 



(4.22) H 2 (A; 0,0) - H 2 (A; 0,1) = i log A - i log A 2 +n 2 d 2 + &rg(A - id) . 

If we use these formulas in (I4.1U|) . JOSJ, l|4.14|l and (|4.15|) we end up with <|2.35[l . 
The asymptotics for f(z) and g(z) are easy to obtain using integration by parts. □ 
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5. Proof of the approximation theorem 

5.1. Main part of proof. We will use the formulas H2.15jl - I|2.17(l for Kg°(u,v). 
Given a sequence {cj}j>i of complex numbers ^ 0, we define the counting function, 

»c(t) = #{j > i; M <t}. 

If n c (t) < Ci 1+l5 for some 5, < § < 1, we can define the convergent canonical 
product 



(5.1) p c (2) = n 1 - 

3=1 V J 

It follows from lemma mi below that £ m , as defined by (|2.43|) is finite. Define 

(5.2) Cj — Cj,77i — Uj+m Urn 5 — ^JjTTO — 2/m Um—j- 

We will first show that 

/ 53 ^ (w + y m )F(w + y m ) = ^ m(M - z) wfiMf t (-t») 

1 ' ' {z + Vm)F{z + V m ) ' zP c {z)P b {-z) ' 

provided z is not a zero of P c or — z a zero of Note that p, and P c are well 
defined by assumption (i) in the theorem. The left hand side of Ij5.3|l is 

(5.4) 

,. w + y m -pr y- -{w + y m Y 
hm —5 t — 

N^oo Z + y m AX y| - (Z + ?/m) 2 

li m w + jfa TT to ~ (j(m + w ))(yj + y m + w) 
n^oc z + y m 11 (y^ - (y m + z))to + y m + z) 

N m-1 N 

lim w + y™ TT ^ -Vm-w _ w yr yj - y m - w -pj- j/j - y m - w 

w^oo z + y m 11^ yj -y m -z n^oo z ,±} N yj - y m - z y j - y m - z 

N+m N—m 

- lim — TT Vm - w TT yj+m _ Vm ~ w 



N~,00 Z AX y m _ y m ,. _|_ ^ XX y j+m _ y?) 



N+m N-m 1 , 

u) TT 1 + TT 1 ~ W / C J 



lim - TT i^i^ TT ±_ 
11 1 + */&,- 11 1 

3=1 



N+m — w( — — ) iV+m „, , „ , N 

to e ^ J TT ^( z -^) Pb(~w)P c (w) 

n^Lz ll e -*(£-£> . 11 * ahow 

j — 1 c j— N— m+l 

which gives the right hand side of l|5.3JI since, Y^jJpp-m+i ;r- ^ as TV ^ oo for a 
fixed m. 

Write v m = v — Cm (we write just m instead of m(a)). Then 

(5.5) |um|<P + A m , |u m |<T + A m . 
Note that 

, , {w + y m - vf (z + y m - uf _ {w- v m ) 2 (z - u m f 

(5 ' 6) 2S 2S " 2S ^—+U(w-z+u-v). 
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If we make the change of variables z 
we obtain, using (|5.3|) . 

(5.7) K%Au,v) 



1 



dw 



z + ijm + u m , w — ► w + y m + v m in (|2.15j) 



a w 2 /2S-z 2 /2S 



(w + v m )P c (w + v m )P b (-w - l> m ) 



z-w + u- w (z + u m )P c (z + u m )P b (-z - u m ) ' 

where we have taken L = a + j/ m + w m . The number a will be specified later and 
satisfies \a\ < X rn . We will choose M = irS/\ m for reasons that will be clear below. 
Note that M > 1 if a is large enough by our assumption on the allowed values of 
S. Below we will need the following estimates of the canonical products. Fix an 
e > 0. There are constants ci, C2 > such that 

(5.8) \P b (w)\, \P c (w)\ < de^^H+H 1 * 5 ), 
for all w £ C and 

(5.9) \P c (z)\ > c-^-^K^'-MA^l^ 
if |Imz| > 1. These estimates are proved using the estimate 

n b {t) < n c (t) < \-H + Ct 1+s , 

see l|5.52[) below, and the following inequalities in 0], p. 19-22. If x = {xk}^ =1 
satisfies n x (t) < Ct 1+S , < 6 < 1, t > and r = \z\, then 

n x (t) 



(5.10) log |P,( Z )|< 8 jr^ r ^ft + r 2 £ 



i 3 



-dt 



for all z G 
(5.11) 



and 



\og\P x (z)\ > -n x (2r)\og(2r) + 



2r 



eft - 8r 2 



2c 



t 3 



dt - 2r 



2r 



S) 



f 2 



dt. 



provided \z — Xj | > 1 for all j > 1. 

Next we will prove an estimate which allows us to restrict the domain of inte- 
gration in l|5.7|l . Fix R > 0. Introduce the following contours: 

(5.12) T a , R : [-R, R]3t^ (j + it, 

T c aR : M\ [-R,R] 3t^ cr + it, 

j MtR : [-R, R}3t->Tt± iM, 

1m,r ■ R \ [~R, R]3t^Tt± iM. 

Let (7,7') denote either (Tct,Jmr) or (r° jj,7Af)- We will show that there is a 
constant cq such that if 



(5.13) 

then 
(5.14) 



( s 

R > Co max < I - — 
V A™ 



1/(1- 



T l+e S 

Am 



1/2" 



5 "' T 



Idzl 



s («, 2 -* 2 )/2S 



1 



\w + v m \ 



P c (w + v m )P b (-w - v m ) 



P c (z + u m )Pb(-z - u m ) 



\z — w + u — v\ \z + u rn \ 
C\ Tt 



< 



5-3/2 
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Write g(t) = X m H + t 1+5 and h(t) = \ m 1 t 1 +^ + t 1+5+e , and Ir = [-R,R]. It follows 
from (|5.8j) and (|5.9jl that the integral in (|5.14jl is 

(5.15) 



<- 



C 
MS 



ds 



dte(- s+ °- t+M V 



■2S 



y/(Tt -<t + u-v) 2 + (ts + M) 2 



x e 



for some constant c > 0. Here we have used \z + u m \ > M and t < exp(</(i)). We 
can now use g(t + s) < 2{g{t) + g(s)) and similarly for h to see that the integral in 
lETTS) is 

(5.16) 



C 

< exp 

" MS y 



ds 



R 1 R 



a 2 + M 2 
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dt- 



cxp 



+ C{g{\a\) + g(\v m \) + h(\u m \) + h{M)) 
(-i±f + C{g{\s\)+h{\t\))) 



y'it-a + u-v) 2 + (s- Mf 



a 2 + M 2 R 2 



- I ds I dl- 

C 



exp i 



+ C(g{\a\) + g{\v m \) + h(\u m \) + h{M)) 
— ^C(g(\s\)+h(\t\)) 



2S 45 

.2 _i_ + 2 



s z +t 



< 



MS 1 / 2 
where 



exp 



- a + u - v) 2 + (s - M) 2 
a 2 + M 2 R 2 



2S 45 

s 2 + 1 2 



+ C(g(\a\) + g(\v m \) + h(\u m \) + h{M)) + A 



+ C(g(\s\VS) + h(\t\VS)) < CS^Xm 1 



A = max 

(s,t)eR : 

since 1 > 5 > X m (essentially). We see that we need 

S " + C(5(H) + 9ilVml) + MM) + KM)) + Csi± ^ X ^ lh ' 

which holds if R satisfies H5.1rt|) . Here we have used \a\ < X m , (|5.5|l and M — 
nS/X m . This proves (|5.14|) . 

We will also need the following estimate. There is a constant C such that 

7r(i/;+u m ) 



(5.17) 



\dw\ 



\dz\ 



7M 



e 2s 



1 



\z — w + u — v\ 



sm 



7r(z+u m ) 



< 



c 



The left hand side of (|5~T7|) is 



(5.18) 

Now, 

(5.19) 



dr / dt 



Tt + u-v + i(a + TTM)\ 



sinY-(v m + g + it) 



sinf^(u m ±i=FM) 



sin - — (v m + a + it) 

X m 



, 2-kt , 22L 
< 2 cosh - — < 4e A " 
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and hence 



Also, by our choice of M, 



sin - — (v m + a + it) 

Am 



< 2e~r< 



sin— (u m ±t^M) 

Am 



(5.20) 

Hence the integral in l|5.18(l is 

sis 

C e 



. , ttM . , ir 2 S 
> sinn — — = sinn 

Am A™ 



(5.21) 



< 



5 sinh 



zdrdt. 



(t - (u- v)) 2 + (t ±<t ± f^) 2 



The contribution to the integral from (i — (it — v)) 2 + (r ± cr ± irS/\ m ) 2 < 1 is 
< C and from the the contribution from the complementary region is < CS. This 
proves H5.17JI . 

It follows from (|5?T4l that 



(5.22) e-^ u -^K* s Ay;u > v) = 



1 



where 
(5.23) 



(27ri) 2 S* 

[w + v m )P c {w + v m )P b (-~w - v m ) 
(z + u m )P c (z + u m )P b (-z - u m ) 

c\ 



dw 



w 2_ z 2 

dze 2S ■ 



7m, r 



z — w + u — V 



s 3 / 



provided R satisfies (|5.13|) . Let aj. = X m k, k > 1. Then 



(5.24) 



zP a (z)P a (-z) = sin — , 

Am. 



and we write 
(5.25) 

(W + V m )Pc(w + V m )Pb{-W - V m ) _ Sm^-(w + V m ) 

(z + u m )P c (z + u m )P b (-z - u m ) sin -£-(z + u m ) 

P c (w + Vm) Pb{-W - V m ) P a {z + U m ) P a {-Z - U m ) 
P a (w + V m ) Pa(—W ~ V m ) P c (z + U m ) P b (- Z - U m ) 

The argument above with Fit) — A" 1 ^ gives 
(5.26) 

1 



sin j^{w + v m ) 
sin y~(z + u m ) 



K *S,i{y} u ' v ) 
where 
(5.27) 



dw 



(2n) 2 S 



dze 2-s 



sin y~ (w + v n 



1M,R 



z — w + u — v sin y~ (z + u m ) 



+n 2 , 



\n 2 \ < ^ e - 2 / 8S 



provided R satisfies (I5.13|) (with an appropriate constant cq that does not depend 
on A m ). Here Kg A {y] u, v) is given by l|2.15|l with L = a + v m . 
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Set 
(5.28) 



n a = 



(2iri) 2 S 



dw 



w 2 -z 2 

dze 2S " 



1 



aia^-(w + v m ) 



7m, r 



z — w + u — v sin y~ (z + u m ) 



P c (w + V m ) Pb{-W - V m ) P a (z + U m ) Pg(-Z - U m ) 
^P a {w + V m ) P a {-w - V m ) P c (z + U m ) Pb(~Z - U m ) 
Then, by E23, (EHH) and i5~2^l 

(5.29) e-^-^K^y; u, v) = K* s ^ u, v)+K 1 -K 2 + K 3 . 

We need an estimate of IZ3. For this we need estimates of 



- 1 



(5.30) 



P c {w + v m ) 



P a (w + v m ) 



1 



and 



P a (z + u m ) 



P c (z + u m ) 



- 1 



and the same expression with b instead of c and a change of sign. We have the 
identity 



(5.31) 
Hence 
(5.32) 



P c {w + V m ) 

P a {w + V m ) 

Pc{w + V m ) 



exp 



< exp 



(w + v m ) 2 n c (t) - n a (t) 



w + v m ~t 
(w + v m ) 2 



(It 



w + v v 



P a (w + v m ) 
Here we can use lemma l5~2*l below with 

(cr + v m f + ,s 2 



t 2 

\n c (t) - n a (t)\ 
t 2 



dt 



(5.33) 



9(t) = 



y/{(7 + V m - t) 2 + S 2 

where w = a+is. Let k m — [v m /X m — 1/2] and choose a = a m = A m (fe m +l/2)— v m . 
Then \a m \ < X m and we have, by (|5.5|) . \a m + v m \ < T + 2X m . By assumption 
T < T Q < 1/4^/77^", and it follows that k m < K = [(2A m% /7y m ")" 1 ] if m (i.e. 

a) is large enough. Using the notation of lemma l5~2l we see that if \t — \ m k\ < o^, 
then 



(5.34) 



t\ > —(\k 



1), 



1 < k < K . Here we have used that ctk < A m /4 if 1 < k < K . If t > K\ m , then 



v m -t\> - 



(5.35) \a n 

if m is large enough. It follows from (|5.34|) and 15.35f) that 

(T + 2A m ) 2 + s 2 



(5.36) 



9{t) < 8 



\ m (\k - k m \ + 1) 



if \t — \ m k\ < afe for 1 < k < K, and 



(5.37) 



<?(*) < 4 



(T + 2A m ) 2 + s 2 



if t > KX m . We ca now use lemma to conclude that 
(5.38) 

1 g(t) K(t) ~ na(t)l dt < [{T + 2A m ) 2 + fl 2 ] f^logC^ 



v rnv 'I in 



-.) + (A m T7„ 



<C[(T + 2A m ) 2 + R 2 ] m — ~ 
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where we have used |s| < R, A m * < Cm° , r] m < Cm 1 + s and A m < C. Hence, by 
P c (w + v m ) 



(5.39) 



- 1 



< exp ( Cm" 11 ^- [(T + 2A m ) 2 + R 2 ] ) - 1 



P a (w + v m ) 

A very similar computation using lemma IS~4l instead gives 



(5.40) 



P b (-w - v m ) 



P a (-w - v m ) 
We also have the estimate 
P a (z + u m ) 



- 1 



(5.41) 
where now 
(5.42) 



P c (z + u, L 



9(t) = 



(z + Um)' 



z + Un 



< exp ( Cm"V- [(T + 2A m ) 2 + R 2 ]) -1. 



f ^ \ n c(t) - n a (t)\ ' 
^ di 



(-.s + u m ) 2 + M 2 
V / (-s + "™-t) 2 + M 2 ' 



< exp 


[/ 




.Jo 



t 



< 



if z belongs to the upper part of Jm,r (the other case is completely analogous). 
We have \s\ < R and 2(R + T + A m ) < \j2- s /r\ m if m is sufficiently large by our 
assumptions on i? and the fact that T <Tq. If t > \j2-Jr\r~n~ we get 



(5.43) 



<?(*) < c 



i? 2 + (T + A m ) 2 
t 



If < t < 1/2-y/Tjm, it could happen that — s + w m is close to t. Here we use a 
similar estimate as above, 



(5.44) 



9(t) < c 



R 2 + (T + \ m ) 2 



with an appropriate fc* (depending on — s + u m ). Using the estimates l|5.43l) and 
(|5.44|l in (|5.41|l and lemma l5~2l we again get 



(5.45) 



P a (z + u m ) 



1 



P c (z + U m ) 

and similarly, using lemma I5~4l instead 

-^~a( ~ Um) 



<exp[Cm-^f^[{T + 2X m ) 2 + R 2 }) - 1, 



(5.46) 



P fc (-z - M m ) 



1 



< exp ( Cm — itt- [(T + 2A m ) 2 + i? 2 



1. 



By our assumptions on i? and T we see that the expression in the exponent in 
(EH! is bounded by a constant. Hence by (|5~i7jl . (|5~2l^l . (|5"59"|) . (|5~IH|l . (|5~l5) and 

we get 



1^3 1 < Cm — r+3-(T 2 +i? 2 



(5.47) 



From (|5~29"|) . (t5~33j> . (|5~77jl and (|5~T7|l it follows that 
(5.48) 



CA ?i 



<- m -^( T 2 + J2 a ) + ^ e -^ 
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Here Kg ^y; u, v) is given by l|2.15[l with L = a + y m + v m and K* s u m , v m ) is 
given by (|2.15() with L = a + v rn . Thus, 

(5.49) 

i p a -\-y m -\-v m -\- M i 



1 

2^1 



J <T~\-ym-\-Vm — M'i 

(T+u m — Mi 



Hence (I5.48f) also gives 
(5.50) 

-U(u- v)KUyjU)V) _ KmUm)Vm) < | TO -^# (T 2 + i?2) + |A™ e -^/85_ 



Note that 



(2-Ki) 2 S J Tl J JM z-w zF(z) 

(2™) 2 S J Vl J 1M z + w zF{z) 

since F(—z) — F(z). We can carry out the same type of computation as above 
to see that (|5.5UI) also holds for K s h {y, u, v). Now Kg h (y, u m , v m ) is approximated 
by Ks(u m ,VM), given by Ks{u,v) = a~ 1 Ls(a~ 1 u,a~ 1 v) and L$ as in H2.20J) . with 
error < -£-, which is smaller than the error term we have in the theorem. This 

— Sa 1 

completes the proof of the approximation theorem. 

5.2. Some lemmas. In the proof above we need some facts about certain numbers 
defined in the theorem. 



Lemma 5.1. The number £ m defined by j2.4°^ is finite. Also if we set £ m — 
Vm — £,m,S, there is, for sufficiently large a, an m = m(a) such that 

(5-51) |Cm(a) - a| < A m ( Q ). 

Proof. In the proof we will need some rather simple facts which we will prove later. 
They are immediate consequences of our assumptions on F. 

(a) F'{t + s) < F'{t) + F'(s), for all t, s > 0. 

(b) F _1 (f + s) < F _1 (t) + F- 1 ^), for all t, s > 0. 

(c) tF'{t) < 4F(t), for all t > 0. 

(d) tF"{t) < F'(t), for all t > 0. 

(e) iF^-^rn) +t) - m- < r] m t 2 for all m > 1, t > 0. 

(e) iFiF-^m) + 1) - m\ < X m H + Ct 1+S for all m > 1, t > 0. 

(f) A m - X m+ j < r) m \ 3 M j, j > 1 

(g) A m+ i < 

Let Cj im and 6 Jim be defined by (|5.2p . Then n c (i) = [F(F _1 (m) + t) — m], where 
[•] denotes the integer part, and hence by (f), 

(5.52) n c (t) < \-H + Ct 1+s 



for t > 0. We have 

oo 



(5.53) ^ = E V( 



1 
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and since F is co nvex, bj, m > Cj, m . Thus < bj^ m -c jtm = 2y m + yj, m -yj+ m < 2y„ 
and we see from ()5.53|) that £ m > and 



(5.54) 



^ = E 



2y n 



7[ (Vj+m - y m ){y m + Vj-m) ' 
Since F(x) < Cx 1+S we have 

(5.55) yj = F- 1 V)>Cj 1 K 1+s > 

and it follows that the series in l|5.54|l is convergent. 

To prove the other statement in the lemma, (|5.51|) . we want to estimate |£„ 
£ m +i|. From (|5.53() we have 



(5.56) f m - Cm+i = E 



Ay m+j - Ay m Aj/ m _ 3 - Ay„ 



where we have used the notation Ayj- — yu+i ~ Uk- If we take t = Ay rn in (e) and 
use (h) we get 

(5.57) |A m -Ay m | < A^ ro 

and together with (g) this gives 



(5.58) 



\Ay m+j - Ay m \ < A m 7? ro (j + 2), 



for j,m > 1. Since A m is decreasing in to, (h) gives 



(5.59) 

Hence, if 1 < j < to, 
(5.60) 



Cj,m Vm+j Vm ^ ,?A 



m+j- 



Here we have used 
A m 

A2m 



(5.61) 



j F'( j F- 1 (2to)) < F'(F- 1 (m) + F~ 1 (to)) <2 



F'(F-!(to)) ~ F'CF-^m)) 
by (a) and (b). Combining l|5.58|l and (|5.6U|) we get 

m— 1 



(5.62) 



E 



Ay m+j - Ay„ 



Cj.m,Cj,m-\-l 



< 4A TO rj m ^ ^-4^ < 36A m ?7 m log to 



if to > 2. From (c), (d) and l|5.55|l we get 
(5.63) 



r ]m =F\y m )<^f±<Cm-^. 



and thus 



(5.64) 



gives 

m — 1 



Ay m+j - Ay m 

i ^7 ,mCj ,m+l 



< C(l0gTO)TO i+^Ar, 
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By (h) and the fact that A m is dereasing we get \Ay m +j — Ay m | < 2A m . Hence, by 



(5.65) 



E 



< 4A r 



^Vm+j - Ay 

7? 



Cj.mCj^rn-i-l 



< 



2A m ^ ^ — — = 2A m / 



t 3 



dt < CX, 



dn c (t) 



Xm.C 



„l-<5 



By (h) and i|5.61[l . c m . m = i/2m — Um > A2 m > ^mA m . It follows that the right 
hand side of lETflSTl is < CA m ((mA^)- 1 + (mAm) 5 " 1 ). Now, by (c), 



F(t) 



> 



and consequently by (|5.55|) 



> Ci 



1-5 



(5.66) 
We find 

(5.67) 



mA,. 



E 



i5 , ( J P _1 (m)) 
^'(F-^m)) 2 

Ay m +j — Ay 



> Cy^" 5 > CmTT* 



l-S 

< Cm 



Cj,mCj',m+l 



Combining (|5.64(l and (|5 . f>T|l we find 

Ay-m+j - Ay 



(5.68) 



E- 



< Cm its A™ . 



J=1 C j,™ c j,™+1 
The same argument that led to 1)5. 64|) gives 

^/2 



(5.69) 



E 



Ay m -j - Ay„ 



< Cm T+^A m . 



Using &j, m &j,m+i > c 3 , m Cj, m+ i and Ay m -j = Ayj- m -i for j > to, we see by similar 
arguments as above that 



(5.70) 



E 



Ay m -j - Ay„ 



j=3m/2 

It remains to consider 

3m/2 

(5.71) £ 

i=m/2 



E 

j=m/2 



Ay„ 



6j,milj,lTi+l 



— 1-15 
< CTO i+^A,-, 



|Ay 



m— j I 



6j,milj,m+l 



m/2 



Ay, 



'm— j,mOm— j,m+l 



Now, by (h), b m _ jjn = y m ~ y i > (to - j)A m and hence, by l|5.66|l. (c) and 
straightforward estimates 



t/2 



(5-72) £ 



A "' <C-^m-^\ rn <Cm~^\ r , 

mXm. 



^m—j,m^m—j,m-\-l 
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Combining this with Ij5.69jl . (|5.70jl and Ij5.71j) we obtain 
Ay m -j - Ay„ 



(5.73) 



E 



^/,m^j,m+l 



< C(logm)m 1+7 A,, 



From 1)5. 56[) . (|5.68|) and (|5.73[l we now get the disired estimate 

(5.74) |£ m - £ m +i| < C{\ogm+l)m-^X m 
for m > 1. Using l|5.74[) and (h) we get, for 1 < r < m — 1, 

m— 1 m— 1 

(5.75) ie m | < & + E - &i - & + c E ( lo s m + i)™ - ^ A ^-i 

3=1 i=i 
< 6 + C'yr + C(logr + l)r"TT7(y ro _i - y r _i). 

Consequently £ m /y m -> as m ^ oo and since 5 < 1, we see that £ m = y m —S^ n 
oo as 772 — > oo . Since Ay m /X m — > 1 as to — * oo it follows from i|5.74[l that 



Cm+l Cm 




Ay m $(£771+1 — Cm) 






Am Am 



as 777, — > oo. Hence |Cm+i — (ml < 3A m /2 if m is sufficiently large. If we take a 
sufficiently large the closest £ m is thus within distance 3A m /4 or 3A m+ i/4, which 
is < X m and we take this m as our m(a). □ 

Our next lemma is 

Lemma 5.2. Let cij = X m j, j > 1 and Cj = Cj, m , j > 1 as above. If g is a bounded 
measurable function on [0,oo) we have 



(5.76) 



g{t) \Mt) ~n a (t)\ ^ 



K r k\ m +a k 



l/\ m Vr, 



\g{t)\dt 



^E 

k=\ 

2 

A m Jl/\ mV „ 



kXm—av. 



t 2 



dt 



\g(t)\ 



dt + C 



dt, 



where a k = X^ n rj m (k + I) 2 , K = [{2X m Jff^) l ] 



Proof. We have n c (t) = [F(t + i f,_1 (m)) - to] and n a (t) = [A~ x i]. The proof of 
(|5.76fl is based on the following claim which we will prove below. 

Claim 5.3. Assume that < t < X m K . 1 < k < K . Then n c (t) — n a (t) if 
\t - X m k\ > a k and \n c (t) - n a (t)\ <lif\t- X m k\ < a k . 



Using the claim we have that the left hand side of l|5.76ll is 



(5.77) 



^E 

k=l 



t 2 



-dt 



2 



^m KQ) -n a (t)\ 
p ti- 



lt follows from (e) that \n c (t) — n a (t)\ < r) m t + 1 and from (f) we get 

|n c (f) - n a (t)\ < n c (t) + n a (t) < ^-t + Ct 1+S . 
Using these estimates in the second integral in (|5.77() we obtain (|5 



□ 
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Proof, of claim HOI If < t < X m K, then t 6 [X m k, X m (k + 1)] for some k, 
< k < K . If furthermore t € [X m k + A m (fc + 1) — a/c+i], then [A _1 i] = k, since 
X^a k = A^?7 m (/c+l) 2 < X 2 m r) m K 2 < 1/4, by our choice of X. Also, A^t-p^t] > 
A" 1 ^ > ry m t 2 since t < X m (k + 1), and A~ x t - ([A~ x t] + 1) < -X^a k +i < ~Vmt 2 
since t < X m (k + 1). Combined with (e) this gives < F(t + F^ 1 (m)) — m — [A~"4] < 
1, i.e. n c (t) = n a (t). On the other hand, if < t < X m K and \t — X m k\ < ctk, then 
n c (t) and n a (t) can differ by at most 1. □ 

We have a similar lemma for n b instead. 
Lemma 5.4. Using the same notation as in lemma \5~^ we have 
(5.78) 



g{t) \Mt)-n a (t)\ dt 



K r k\ m +u k 



+ r n m4 + + f r m + c r m*, 

Jl/4^ r) m/2 1 A m Jm\ m /2 1 Jm\ m /2 1 

where the second integral in the right hand side is present only ifmX m /2 > 1 /4^/7y TO / 2 . 

Proof. An integration by parts and the fact that F" is decreasing gives 

\m - F(F- l (m) - t) - X m H\ < 2rj m/2 t 2 

if < t < y rn - y m /2- Note that by (h), y m - y m/2 > mX rn /2. Now, n b (t) equals 
[m — F(F~ 1 (m) — t)] if < t < y m and m + — y m ) if t > ?/ m . It is clear that 
nb(t) < X~^t for < t < y m . If t > y m , then using < F(t) we find 

n 6 (f) < + F(t - F _1 (m)) < F(t) < Ci 1+l5 . 

The proof of (|5.78|) now proceeds in the same way as the proof of lemma □ 

It remains to prove the statements in the begining of this subsection, (a) and 
(b) are immediate consequences of our assumptions on F. To prove (c) write 

F(t)= f F'{s)ds> f F'(s)ds> t -F'(i), 

<J J t j 2 

since F' is increasing. By (a), F'(t) < 2F'(t/2), and (c) follows. The statement (d) 
follows from the fact that F" is decreasing. To prove (e) write 



(5.79) iFiF-^m) +t) -m- X m H\ 



■ F^ 1 (m)+t 

{F- l {m) + t- s)F"(s)ds 
< F"{F~ l (m))t 2 = r] m t 2 , 



since F" is decreasing. For (f) write 

(5.80) F(F~ 1 {m)+t)-m= I F' '(s)ds < tF 1 '(F -1 (m) + t) 

J F- 1 (m) 

< tF'dF-^m)) + tF'(t) < X m H + Ct 1+5 7 
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by (a), (c) and the fact that F' is increasing. To prove (g) write 

(5 81) A A ■ — F'(F-H™ + j))-F'{F-Hm)) < 1 f Vm+i F „ (s)ds 
(5.81) A m - X m+J FI{F - 1{m))F , {F - 1{m + j)) < F 

S.2 /„, . . „, \_„ \2 / °± - \3 



Finally, to prove (h) we write 

r m+1 dt 



f 

Vm+i -y m = 

J m 



Since F' o F^ 1 is increasing the right hand side is < A m and > A m+ i. 
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